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Abstract 

The contribution of virtual excitations to the energy levels of electronic and muonic hy- 
drogen is investigated combining a model-independent approach for the main part with 
quark model predictions for the remaining corrections. Precise values for the polarization 
shifts are obtained in the long-wavelength dipole approximation by numerically integrat- 
ing over measured total photoabsorption cross sections. These unretarded results are 
considerably reduced by including retardation effects in an approximate way since the 
average momentum transfer (together with the mean excitation energy) turns out to be 
larger than usually assumed. Transverse and seagull contributions are estimated in a sim- 
ple harmonic oscillator quark model and found to be non-negligible. Possible uncertainties 
and improvements of the final results are discussed. 
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1. There has been tremendous progress in the laser spectroscopy of hydrogen and deuterium 
atoms which now are even sensitive to small nuclear and proton structure effects. One 
of these - traditionally the least understood - is the virtual excitation of the nucleus which 
acts back on the bound lepton. While the effect in deuterium is comparatively large and 
has been evaluated theoretically with increasing sophistication and reliability 0, the proton 
polarization shifts have not received much attention up to now. Khriplovich and Sen'kov || 
have estimated the shift in the electronic 15-state as —71 ± 11 ± 7 Hz using values for the 
static proton polarizabilities and assuming a mean excitation energy of ~ 300 MeV. They 
attribute the quoted errors to the use of a relativistic approximation for the electron and to 
the experimental values of the polarizabilities. However, experience gained in many decades of 
nuclear polarization calculations has told us that the use of an average excitation energy (the 
"closure approximation") can be a considerable source of uncertainty unless one calculates it 
precisely. In addition, one has to make sure that other ingredients to the polarization shift 
(higher multipoles, transverse excitations) are well under control before a definite answer can be 
given. It is the purpose of this note to re-evaluate the polarization shift without the questionable 
use of a mean excitation energy and other simplifying assumptions. Since an experiment is in 
progress at PSI to measure the Lamb shift in muonic hydrogen |4[] we will also evaluate the 
polarization shift for this case. Actually, with the experimental accuracies achievable in the 
near future, it turns out that the proton polarization shifts are of much greater importance for 
muonic than for electronic hydrogen. 

2. For light electronic and muonic atoms the energy shift due to virtual excitations can be 
written as an integral over the forward virtual Compton amplitude. This quantity in turn 
may be expressed by its imaginary part, i.e. the structure functions Wx^i^^Q 2 ) which are 
measurable in the inclusive reaction e + p — > e' + X. In the absence of detailed experimental 
information for all relevant values of momentum transfer Q and energy transfer v it is customary 
to apply the long wavelength (or unretarded dipole) approximation which should be valid for 
Q (r 2 } 1/2 <C 1 where (r 2 ) 1 ^ 2 is the root-mean-square radius of the proton and Q an average 
momentum transfer. In this limit it is possible to express the structure functions by the 
experimentally measured photoabsorption cross section <t 7 (z/). Bernabeu and Ericson have 
derived in this way the following expression for the energy shift || 
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Here a is the fine-structure constant, m the lepton mass and ip n i(0) the lepton wave function 
at the origin which is only non-zero for 5-states. The correction factor accounting for the 
variation of the wave function over the proton radius can be safely neglected at this level of 
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accuracy. We first note that it is possible to perform all integrations in eq. (0) [] and to obtain 
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The function f(x) has quite different limits for small and large arguments: 
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which makes a crucial difference between nuclear polarization shifts, say, in muonic deuterium 
and in electronic hydrogen. Krhriplovich and Sen'kov have only retained the leading logarithm 
in eq. (|6|) which is justified considering the much greater error which comes from pulling out 
the logarithmic term evaluated at the excitation energy of the A (1232) and expressing the 
remaining integral in terms of the sum of polarizabilities 
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3. In view of the well-known deficiencies of the closure approximation we have decided to 
evaluate the polarization shift by integrating numerically over the experimentally measured 
photon absorption cross section. We have taken the recent Mainz data from v = 200 — 800 
MeV, the older Daresbury data || from v = 800 — 4215 MeV and the parametrization 
a 7 (v) = 96.6 /ib + 70.2 fib GeV 1/2 /^ above 4215 MeV. Below 200 MeV four angular dis- 

) have been converted to total cross sections 
) = 18.3/ibMeV" 1/2 yjv - u th . The 
numerical integration was done by simple Simpson integration over linearly interpolated data 
points and the whole procedure was checked by evaluating eq. (0). We obtained 



tributions in 7p — > 7i + n (from Fig. 1 in ref. 
and the threshold behaviour has been parametrized as a. 
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in good agreement with a recent analysis JTIJ but smaller than the value 14.2 • 10 4 fm 3 used 
in Ref. M. This is due to the fact that the older cross sections which have been analysed in 



Ref. |yj are systematically higher than the new Mainz data. 

The unretarded results for the electronic (e) and muonic (//) polarization shifts are then 
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This has also been noted before 
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In view of the small electron mass it is not surprising that the relativistic approximation (^j) 
(including the constant term) agrees with the exact result to more than 4 digits. Equation 
(0) is larger than the estimate of ref. || because the integral over virtual excitations gets 
contributions well above the A-resonance. This can be easily seen from any graph of the total 
photoabsorption cross section versus photon energy but made more quantitative by asking for 
the value of the mean excitation energy v which - when substituted into the logarithm of 
ec L- @ - gives the same result for the shift. We find v ~ 410 MeV both for electronic and 
muonic hydrogen. Therefore predictions assuming that only the A isobar contributes to virtual 
excitations of the proton |12|] are not very reliable. In the muonic case the approximation (|6|) still 



gives more than 96 % of the exact numerical result whereas the non-relativistic approximation 
(0) overestimates it by more than 40 %. This is, of course, consistent with the fact that the 
mean excitation energy is nearly four times larger than the muon mass. 

4. The large value of the mean excitation energy also casts some doubt on the use of the 
unretarded dipole approximation. For example, in a simple constituent harmonic oscillator 



quark model one naively expects an associated mean momentum transfer of Q ~ yM quar kZ/ 
350 MeV. This will have an appreciable effect when inserted in the elastic form factor which, e. g. 
in the dipole approximation is given byF (Q 2 ) ~ 1 /(1+Q 2 /Q 2 ) 2 with Ql = 0.71GeV 2 . Inelastic 
transition form factors to low-lying resonances have approximately the same Q-dependence 
apart from threshold factors which are characteristic for the specific angular momentum of the 
resonance (see below for a non-relativistic example). Therefore one obtains a rough estimate 
for the effect of retardation when the square of the elastic form factor Fq(Q 2 ) is inserted in 
the Q 2 -integral of eq. (0). Instead of the dipole form factor we have employed a more realistic 
parametrization of the charge form factor of the proton given in eq. (9) and table 3 of ref. 



I"3f . A careful numerical evaluation of the remaining double integral then gives for the retarded 



polarization shifts 

88 9 

*Ks = ~ Hz (11) 
= (12) 

Use of the dipole form factor changes the numerical values to 89.3 and 114 for electron and 
muon, respectively. Equation ([TT|) is nearly 20% smaller in magnitude than the unretarded 
result and the polarization shift in muonic hydrogen gets almost halved due to retardation 
effects. This reduction can be translated into a mean momentum transfer which - when 
inserted into the square of the elastic form factor - cuts the unretarded values by just this 
amount. In this way one obtains Q ~ 180 and 340 MeV/c for the electronic and muonic case 
respectively. 

5. At present this seems to be the best mo del- independent estimate for the polarization shifts in 
hydrogen but it still neglects higher multipoles and relies on a standard but not very well tested 
procedure to correct the unretarded dipole approximation. Contrary to the deuteron case where 
excellent potential models are available, a reliable (relativistic ) model of the nucleon does not 
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exist to calculate all these contributions. Here we take the simple non-relativistic harmonic 
oscillator quark model to estimate them. In this model it is easy to evaluate analytically the 
longitudinal (inelastic) structure function 

^■«>%| < i3 > 

as well as the transverse one 
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Here y = (qb) 2 /3 where b is the oscillator length and uj = l/(M quar k& 2 ) the harmonic oscillator 
frequency. In this model the excitation spectrum consists of sharp lines at Nu shifted by 
the recoil energy of the proton and the elastic form factor is gaussian (F^q) = exp (—y)) so 
that the proton rms radius is directly given by the oscillator length b. It should be noticed 
that the first and the second term in eq. fllip come from excitations by the spin and the 
convection current, respectively. Note also that the above structure functions fulfill Siegert's 
theorem lim g ^ Sl(v, q)/q 2 = lim g ^ St(v, <?)/(2^ 2 ) = o ^{y) j (Am 2 av) and that all excitations 
are indeed multiplied by the square of the elastic form factor. In the low-g limit the first excited 
state exhausts the dipole absorption cross section which is not a very realistic feature of the 
model. Other shortcomings are the well-known inability |L5| of the harmonic oscillator quark 
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model to reproduce both the empirical rms-radius ( (r 2 ) = 0.86 fm ) and the polarizabilities 
(0) if the constituent quark mass is fixed to M quar k = M prot on/3. We also make this choice 
because the masses of the constituents should add up to the total mass in a consistent non- 
relativistic treatment. In addition, this value gives the correct recoil energy and leads to 
a reasonable result for the magnetic moment of the proton. This leaves only the harmonic 
oscillator length as free parameter which we have fixed in such a way that the results from 
the unretarded dipole approximation are obtained. In this way the quark model results can 
be directly compared with the shifts evaluated in the model-independent approach. It is quite 
obvious that a non-relativistic model becomes inadequate for excitation energies and momenta 
of the order of the constituent mass. However, here our aim only is to obtain a rough estimate 
for the remaining corrections beyond the retarded dipole approximation and for this purpose 
the harmonic oscillator quark model may be not totally useless. 

To obtain a quantitative estimate we have inserted the analytic expressions ( |T3"| , O ) into 
the formulae in the Appendix of ref. |I3 (which keep the relativistic kinematics for the lepton), 
summed over oscillator shells up to iV max = 20 and integrated numerically over the momentum 
transfer up to g max = 3000 MeV. Special care has to be exercised because of the apparent 
singularity in the transverse weight function at q = which is canceled by the seagull contri- 
bution. The latter one is required for a gauge-invariant treatment of non-relativistic systems. 

2 The structure functions S^/t are linear combinations of the usual W1/2 and are more convenient in the 
framework of ref. [[u] which uses the energy transfer v and the three-momentum transfer q = |q| as variables 
whereas ref. H employs the invariants v and Q 2 . 
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Other numerical difficulties arise from the very different scales of electron mass and mean ex- 
citation energy which requires very high numerical accuracy (up to 60 x 72 gaussian points 
compared to 6 x 72 in the muonic case) and from the slow convergence (~ 1/-/V^ ax ) of the sum 
over excitations for the spin current which peaks at high momentum, i.e. high N . The latter 
problem has been overcome by an analytic resummation and the overall numerical stability 
has been checked by varying g max , N max and the number of gaussian integration points. The 
results of the calculation are collected in Table 1. As expected the b- values are too small to 
account for the proton radius while the upvalues (~ 280 MeV) seem reasonable for low-lying 
nucleonic excitations. A simple way to cure these deficiencies is to attribute the gaussian form 
factor to the quark core and to introduce an extra formfactor 1/(1 + f3 2 q 2 ) for the meson cloud 
which surrounds the core and brings the rms-radius of the model in agreement with the ex- 
perimental value nj. Although this is rather ad hoc and theoretically not very appealing we 
have included this variant also in Table 1. Note that the unretarded dipole approximation 
(and therefore the value obtained in ref. ||) also includes some transverse excitations : eq. 
(§) would read (In Ax + l)/x 2 if only longitudinal excitations are kept. While the longitudinal 
excitations (including all multipolarities) are seen to dominate, virtual transverse excitations 
induced by the currents cannot be neglected since the mean momentum transfer is of the order 
of the constituent quark mass. This is particularly important for the spin current in muonic 
hydrogen because its contribution grows with momentum transfer. 





electronic hydrogen [Hz] 
(A) (B) 


muonic hydrogen [fieV] 
(A) (B) 


unretarded dipole (input) 
retarded dipole 


- 106.4 
- 94.3 - 92.1 


- 201.8 
- 131.6 - 121.1 


full longitudinal 
spin current 

convection current + seagull 
total 


- 78.2 - 75.9 
- 7.2 - 3.8 

- 19.2 - 18.9 
- 104.6 - 98.6 


- 120.2 - 108.8 

- 39.1 - 20.9 

- 16.9 - 15.4 

- 176.2 - 145.1 


correction to retarded dipole 


- 10.3 - 6.5 


- 44.6 - 24.0 



Table 1 : Polarization shifts to the IS" level in electronic and muonic hydrogen evaluated in 
the harmonic oscillator quark model. The parameters (M quar k = 312.8 MeV, b = 0.657 fm for 
electronic hydrogen, b = 0.674 fm for muonic hydrogen) have been fitted to give the unretarded 
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dipole approximation and the corresponding results are given under the heading (A). In case 
(B) an additional meson-cloud form factor has been introduced to reproduce the experimental 
proton radius. 

The ad-hoc introduction of the meson-cloud form factor reduces all contributions and brings 
the retarded dipole approximation more in accord with the calculation employing realistic form 
factors. 



6. In conclusion, we have evaluated the proton polarization shifts in electronic and muonic 
hydrogen in a fairly model-independent way by integrating over the experimental photoabsorp- 
tion cross section and accounting for retardation by use of the empirical elastic form factor. 
The remaining contributions (mostly from transverse excitations) have been estimated in a sim- 
ple harmonic oscillator quark model and are therefore rather model-dependent and uncertain. 
Since it is physics in the resonance region which dominates these contributions, the theoretical 
situation will probably remain so unless better experimental information from inclusive ep - 
scattering in this region is available. For our final values we add the corrections (B) listed 
in Table 1 to eqs. (|TT| , |1^ ) and assign an error to them which covers the values obtained in 
case (A). This seems reasonable and prudent in view of the mentioned uncertainties and the 
inadequacy of the harmonic oscillator quark model. In addition, we take the difference between 
the retarded dipole result obtained with realistic form factors and the one with the gaussian 
form factor for the core and monopole form factor for the meson cloud as error estimate for 
the model-independent contribution. Adding the errors linearly we obtain in this way our final 
result for the polarization shifts 

qc; 4- 7 

*Ks = ^ Hz ( 15 ) 

. n 136 ±30 , r 

A<s = ^— /ieV . (16) 

The first value is one order of magnitude below the present experimental accuracy (840 Hz) in 
the 2S — IS transition H whereas the planned Lamb shift experiment in muonic hydrogen @ 
aims for a precision which is comparable to the uncertainty in AE^g. Incidentally, the proton 
polarization contribution to this Lamb shift has nearly the same magnitude as the hadronic 



vacuum-polarization correction [[L7| which, however, is more precisely known. A more accurate 
evaluation of the former contribution is therefore needed for a better determination of the 
proton radius from the muonic Lamb shift experiment. 

Note added: After submission of the manuscript additional calculations of the muonic polar- 
ization shift using different methods have been reported. Faustov and Martynenko ]TB| obtain 
nearly the same value as reported in this paper whereas Pachucki's number [JOI] is slightly 
lower. 
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